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Abstract 

I We analyze the cosmology of the Randall-Sundrum model and that of compact 

i ' brane models in general in the presence of a radius stabilization mechanism. We 

find that the expansion of our universe is generically in agreement with the expected 
effective four dimensional description. The constraint (which is responsible for the 
appearance of non-conventional cosmologies in these models) that must be imposed on 
the matter densities on the two branes in the theory without a stabilized radius is a 
consequence of requiring a static solution even in the absence of stabilization. Such 
constraints disappear in the presence of a stablizing potential, and the ordinary FRW 
(Friedmann-Robertson- Walker) equations are reproduced, with the expansion driven 
by the sum of the physical values of the energy densities on the two branes and in the 
bulk. For the case of the Randall-Sundrum model we examine the kinematics of the 
radion field, and find that corrections to the standard FRW equations are small for 
temperatures below the weak scale. We find that the radion field has renormalizable 
and unsuppressed couplings to Standard Model particles after electroweak symmetry 
breaking. These couplings may have important implications for collider searches. We 
comment on the possibility that matter off the TeV brane could serve as a dark matter 
candidate. 
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1 Introduction 



It has been widely understood over the past two years that a very promising route towards 
reconciling the apparent mismatch of the fundamental scales of particle physics and gravity 
is to change the short distance behavior of gravity at distances much larger than the Planck 
length |], H, H|]. One prominent suggestion has been to lower the fundamental scale of gravity 
all the way to the weak scale by introducing large extra dimensions PJ. This possibility 
attracted much attention; in particular, the presence of large extra dimensions opens up the 
possibility for new cosmological scenarios for the early Universe, which has been the subject 
of study of several interesting papers [4-11]. While the possibility of the existence of large 
extra dimensions is very exciting on its own, the existence of a big hierarchy between the 
weak and the Planck scales requires the radius of the extra dimension to be much larger 
than its natural value (see [0 for a scenario in which this might naturally occur). 

Randall and Sundrum (RS) suggested a different setup [13, 14], in which the extra dimen- 
sions are small, but the background metric is not flat along the extra coordinate; rather it is 
a slice of anti-de Sitter (AdS 5 ) space due to a negative bulk cosmological constant balanced 
by the tensions on the two branes of this scenario. In this case, the curved nature of the 
space-time causes the physical scales on the two branes to be different, and exponentially 
suppressed on the negative tension brane. This exponential suppression can then naturally 
explain why the physical scales observed are so much smaller than the Planck scale. For 



early papers related to this subject see |T5], |i6[ ; generalizations of the RS models can be 



found in [17-26], embeddings into supergravity or string theory are discussed in [27-32]. 
Some of the aspects of the cosmology of these models have been examined in [33-39], and 



phenomenological consequences in |40j, [41 



Clearly, the cosmology of this model can be very different from ordinary inflationary 
cosmology in four dimensions. A more detailed study of the early cosmology of the Randall- 
Sundrum model and that of brane models in general (including the large extra dimension 
scenario of JIJ) has been however hindered by one obstacle: it seemed that the late cosmology 
of brane models will generically deviate from the usual FRW Universe in the 4D theory on 
our brane. This would bring into question most of the results of the early papers about the 
cosmology of brane models, where it has been commonly assumed that the late cosmology 
(from BBN to present) of these models is given by the standard FRW cosmology. This 
conclusion is reached by applying the results of f33fl , who examined the solutions to Einstein's 
equations in five dimensions on an S 1 /Z 2 orbifold, with matter included on the two branes, 
and with no cosmological constants in the bulk or on the brane. The result of this study was 
that the expansion of the brane does not reproduce the conventional FRW equations. This 
result has been applied to the RS model in |35], [3(| , and it was found that the cosmology of the 



model with an infinitely large extra dimension (where one lives on the positive tension brane, 
the "Planck brane" ) in fact reproduces the ordinary FRW Universe. However, for the model 
which solves the hierarchy problem (in which our universe lives on the negative tension brane, 
the "TeV brane" ) there is a crucial sign difference in the Friedmann equation. However, this 
conclusion about the wrong-signed Friedmann equation (and the general conclusion that 
brane models have unconventional late cosmologies) has been reached in a theory without 
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an explicit mechanism for stabilizing the radius of the extra dimension, the radion field. It 
was already pointed out in |19|, |33], |35| that the effect of a radion stabilizing potential could 
significantly alter this picture. 

In this paper we confirm this intuition. We examine the cosmology of brane models, in 
particular that of the RS model after the radius is stabilized (for example by the mechanism 
suggested by Goldberger and Wise [|HJ or Luty and Sundrum p0[] , but the details of the 



stabilization mechanism will be irrelevant). We explain that the origin of the unconventional 
cosmologies is not due to a breakdown of the effective 4D theory (neither is it due to the 
appearance of a negative tension brane in the case of the RS model) , but rather to a constraint 
that the matter on the hidden and visible branes are forced to obey in order to ensure a static 
radion modulus. This constraint ensures a static solution for the radius even in the absence 
of a stabilization mechanism. Once such a constraint on the matter content appears in a 
theory, all the results about non-conventional cosmologies mentioned above follow. However, 
if the radius is stabilized, such constraints disappear (they turn into an equation determining 
how the radius is shifted from its minimum due to the presence of matter on the branes), 
and the ordinary FRW equations are recovered for the effective 4D theory for temperatures 
below the weak scale. We also examine the kinematics and dynamics of the radion field of 
the RS model in the presence of a stabilization force. We find that its natural mass scale 
is of the order of the weak scale, and that the shift in the radius is tiny for temperatures 
below the TeV scale. We also find that the radion TeV suppressed interactions with the 
Standard Model (SM), which results in renormalizable couplings to Standard Model fields 
after electroweak symmetry breaking. The existence of these couplings ensures that the 
radion can decay before nucleosynthesis, and may also have important consequences for 
collider searches. 

We also comment on the possibility of matter in the bulk or on the Planck brane serving 
as dark matter candidates. Matter in the bulk naturally has TeV suppressed couplings to SM 
fields, so its annihilation cross section may result in interesting relic densities today. Matter 
on the Planck brane was presumably never in thermal equilibrium after inflation, and it is 
therefore not clear whether it could have interesting relic densities. 

This paper is organized as follows: in Section 2, we present the effective 4D description 
of the RS model. We explain the physical reasons why the constraint (which we derive 
in Appendix A together with an approximate solution in the bulk for the case without a 
stabilized radion) should become irrelevant in the presence of a stable radion and find the 
effective 4D action for the radion-graviton system. In Section 3, we present a detailed solution 
to the 5 dimensional Einstein equations confirming the physical intuition of Section 2: once 
the radius is stabilized, the ordinary FRW Universe is reproduced. A similar solution to the 
case with vanishing background cosmological constants is derived in Appendix A, illustrating 
that our arguments are generic for brane models. In Section 4, we discuss the kinematics 
and dynamics of the radion field based on the effective 4 dimensional action obtained in 
Section 2. We also find the couplings of the radion field to the Standard Model particles. 
In Section 5, we discuss the issues related to (dark) matter on the hidden brane or in the 
bulk. We conclude in Section 6. Appendix A gives an approximate solution to the Einstein 
equations for the RS model with matter on the branes, but no stabilization mechanism for the 
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radion. Appendix B discusses the cosmology of brane models with vanishing cosmological 
constants. In Appendix C, we derive the radion mass in the Goldberger-Wise model for 
radius stabilization. In Appendix D, we show that the Hubble law in the effective theory is 
consistent with the Newtonian force law between particles on the TeV and Planck branes, 
while in Appendix E we calculate the interaction strength of matter in the bulk to matter 
on the Planck and the TeV branes. 



2 The Effective Four Dimensional Theory 

In this section we construct the effective 4 dimensional equations of motion of the RS model. 
We demonstrate, that the constraint on the matter on the two branes (which is derived 
explicitly in Appendix A) in the theory without a stabilized radius can also be obtained 
from the effective theory. This constraint is a consequence of requiring a static solution to 
the equations of motion without a stabilization mechanism. 

Throughout this paper (except in Appendix B) we consider the RS model perturbed by 
matter on the two branes. The metric is given by 

ds 2 = n(y, tfdt 2 - a(y, t) 2 (dx{ + dx\ + dx\) - b(y, t) 2 dy 2 , 

= g AB (x,y)dx A dx B . (2.1) 

The two branes are located at y = and at y — 1/2. The Einstein tensor for this metric is 
given by 




Here primes (dots) denote derivatives with respect to y (t). The Einstein equation is given 
by G a b = K 2 T ah , where T ab is the energy-momentum tensor, and k 2 = where M is 

the five dimensional Planck scale. There is a contribution to T from the bulk cosmological 
constant of the form 

TT = 9abK (2.3) 
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and from the branes 

rpb,brane = M diag + p ^ ^ - p„ V, - p„ K - P*, 0) + 

^~^ diag (- V + p, -V - p, -V -p,-V- p, 0), (2.4) 
b 

where K is the (positive) tension of the "Planck" brane at y — 0, p* and are the density 
and pressure of the matter included on the positive tension brane (where we assume an 
equation of state = w*p*) as measured with respect to the metric g, and p and p are the 
density and pressure of the matter on the negative tension brane (the "TeV" brane), again 
measured with respect to g. Hereafter quantities measured with respect to the metric g are 
refered to as "bare " quantities. We know that in the limit p,p, p* , p* — > we want to recover 
the static Randall-Sundrum solution^ of the form 

n ( y ) = a {y) = e - |y|m ° 6 °, (2.5) 

where the relations between A, V*, V and are given by 

K z 

A = (2.6) 
The effective 4D Planck scale is then given by 

(SvrG.v)- 1 = Mp t = i^?, (2.7) 

where the notation for the present-day value of the warp factor 

Q Q = e - mobo/2 (2.8) 

has been introduced. 

Investigation of the cosmology of brane models has shown that in order to find solutions 
to Einstein's equations with matter on the branes there appears to be a constraint between 



the matter on the two branes, if one requires that the extra dimension remains static [|33], [34 
The appearance of such constraints, as explained in Appendix A for the case of the RS model, 
generically leads to non- conventional cosmologies. The constraint between the "bare" matter 
density p* on the Planck brane, and "bare" matter density p on the TeV brane for the case 
of the RS model is found in Appendix A to beQ 

p* = - P n 2 . (2.9) 



*Note that in this paper we use a different notation than Refs. jl3|, [TJ]. The parameter k of (T^, |l4| is 
denoted here by mo, while the radius of the extra dimension r c is denoted by bo- For the coordinate along 
the extra dimension we use y instead of and for the 5D Planck scale we use M 3 instead of 2M 3 . In our 
notation k 2 = 1/2M 3 . We also place the branes at y = and at y = ^ instead of and w. 

TFor the case of one extra dimension with two branes at y = and y = 1/2 and with vanishing cosmological 
constants this constraint is of the form /?*a(0) = — pa(h) p3|. 
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In addition, the computation in Appendix A of the Hubble parameter, H, for the induced 
metric at the Planck brane implies that p* > 0. So the above constraint implies that the 
energy density on the TeV brane must be negative. Since this is at odds with phenomenology, 
it is important to understand the physical origin of these constraints before concluding that 
this model (and brane models in general) is phenomenologically unacceptable. 

In Appendix A it is shown that although there is a constraint between the Planck brane 
and TeV-brane energy densities, an effective 4D theory does exist. For example, for small 
perturbations in the energy densities the two branes expand at practically the same rate, 
and by an amount that agrees with the effective theory expectation. This is consistent with 
what one expects in a 4D effective theory where there is a (approximate) uniform expansion 
of the two branes. This suggests that the origin of the unconventional cosmologies is the 
above-mentioned constraint, rather than the fact that the TeV brane has negative tension, 
or that there is a sick 4D effective theory. 

In this section the 4D effective theory is constructed from two equivalent approaches. 
First, we directly average the 5D Einstein equations over the bulk to obtain some 4D equa- 
tions. Alternatively, we average the 5D Einstein action over the bulk to obtain a 4D effective 
action. From both approaches we obtain the following picture: without a radion potential 
and for generic energy densities on the two branes, the radion runs off to infinity. This can be 
avoided, again without a radion potential, by tuning the two energy densities in precisely the 
manner prescribed by Eq. (|2.9| ). Thus the above constraint equation is just a consequence 
of requiring that the radion modulus remains static in the absence of a stabilizing potential. 
Once the radion is stabilized the above constraint becomes irrelevant, and one obtains the 
conventional late cosmology. 

2.1 The Averaged Einstein Equations 

In this subsection we average the 5D Einstein equations over the bulk in order to demonstrate 
that without a stabilizing potential the system is over-constrained once we require that the 
radion modulus is static. To perform the averaging, we linearize the metric about the RS 
solution: 

a(y,t) = a{t)Q{y,b(t)){l + 5a(y,t)) 
n(y,t) = n(y,b(t))(l + 5n(t,y)) 

b(t,y) = b(t)(l + Sb(y,t)). (2.10) 

Here we keep the time dependence b{t) in the warp factor Q, and we will use the notation 

tt = tt(y, b(t)) = e - m ° myl/2 , n b = fi(l/2, b(t)) . (2.11) 

The value of Qb, when b = bo = constant, is then given by Qq. In this expansion we 
allow 5fi = 5a, 5n, 5b to be of the form 5f = 5 f (p(t) , p*(t) , y) and we only assume that 
5f ~ 0(p, p*). This last assumption is reasonable, for in the limit that p, p* — > we should 
recover the RS solution. Thus an expansion in 5fi is equivalent to an expansion in p, p*. It 
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is then sufficient to work to linear order in Sfi. Terms such as dSa/dt ~ p ~ p 3 / 2 are higher 
order in p, p* and are neglected, while we work to all orders in b(t) and its time derivatives. 
We also include a radion potential 

a 3 V r {b) = -b{t) J dy^CR = b[t) J dyO 4 T ° (2.12) 

in the computations presented here, and set V r = when desired. Here Cr is some non- 
specified bulk dynamics responsible for generating a potential. Although our focus in this 
section is the constraint obtained without a radion potential, we include it here for later 
application. 

Recall that in classical electromagnetism on a manifold without boundary one integrates 
V • E = p over the manifold to conclude that the sum of the charges must vanish. So here, 
to see if there is a topological constraint on the component energy densities, we might try 
integrating the analog of Gauss' Law in Einstein's theory. Namely, consider 

/ dy fi 4 G° = k 2 J dy fi 4 T ° . (2.13) 

We shall see that this equation does not lead to a topological constraint. Rather, this 
equation, combined with the average of the other Einstein equations demonstrates that the 
constraint on the energy densities follows from requiring a static extra dimension even when 
there is no radion potential. 

Substituting the above expansion, Eq.( |2.10| ), into Eq. ( |2.13| ) and integrating gives 



? + '^T^f Si " ^T^f h = ^T^f (* + <« + W» + * (2.14) 

where e 2 = 0((5a) 2 , (8h) 2 , {5b) 2 ). Note that there are no corrections to H 2 linear in the 
perturbations. Thus the corrections to the Hubble formula from the details of the bulk 
geometry are important only at high temperatures. Also note that this equation reduces to 
the conventional FRW Hubble law when the energy density in the radion is small (i.e. b 
and V r are negligible) compared to the other sources. (Recall that 87iGn = K 2 m /(1 — ^o))- 
Further, in this limit the expansion rate is set by p* + pfi^, which is the sum of the physically 
measured energy densities on the Planck and TeV branes, respectively. For small oscillations 
of the radion, the expansion of the universe is ordinary FRW. Thus the bulk average of the 
G o equation resulted in the 4D Hubble Law. 

Repeating the above averaging procedure for the Gij equation gives, for 6 = : 

• 2 

°_ + 2 - = -SvrGV (p. +p0 4 - V r (b )) + t 2 (2.15) 
a z a v J 

This is just the FRW "pressure " equation. As with the G 00 equation, there are no corrections 
of 0(e). 

The unaveraged linearized G55 equation is given by (again for b = 0) : 

3 ^ (3( W + 5n >) _ 3 J. (X + ^ = -n 2 b 2 (2Mb - Tf) + e 2 . (2.16) 
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Here we do not average over the fifth coordinate. If this is performed, then on the LHS 
quantities such as 5a and 5n appear, whose value depends on the detailed form of the solution 
in the bulk. The results of Appendix A could then be used to obtain the same constraint 
obtained below in Eq. (|2.21| ). As seen below, more generality is obtained by not averaging 
the G55 equation. Next the "jump equations" for 5a' and <W, derived in |33[, given by 

6a% = -ln 2 b(p, + 5bV,) ,6a'\ 1/ 2 = l^b(p + 5bV) , (2.17) 
o o 



5n'\ = -K 2 b(3p, + 2p, -5bV*) , 5n'\ 1/2 = --n 2 b (3p + 2p - 5bV) , (2.18) 
o o 



are inserted into Eq. fl2.16|) to obtain (for example at y = 0) 



1 K 2 m (-p* + 3p*) - p- + - j = -^ 2 T 5 5 | . (2.19) 
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In the following we show that without a radion potential the three equations (|2.14j ), ( |2.15|) 
and ( [2.19 ) imply a constraint between the matter densities if a static solution (b = constant) 



is imposed. The important point about Eq. ( |2.19| ) is that in the absence of a radion potential 



the system is over-constrained. To see this, we can eliminate the a(t) dependence by using 
the other two equations, Eqs. Q2.14| ) and ( |2.15| ). For a static solution for b, i.e. 5b = and 



bo = 0, and with V r = T 55 = 0, this gives 

- -U 2 m (-p* + 3p*) = (p» - + (p - 3p)Qq) (2.20) 

O O 1 — ' 

which simplifies to 

(-3p + p)n* = (3p*- P *)n 2 . (2.21) 

In order that this equation remains consistent with the two conservation of energy equations, 
Eqs. ( |A.7|) and (|A.8|) , a further fine tuning (for w 7^ 1/3) w = is needed. So the above 



constraint is then 

p* = -Sl 2 p , (2-22) 

which is the same constraint obtained from the explicit solutions Appendix A, Eq. (|A.13|) . 
Alternatively, the same constraint could have been obtained by comparing the unaveraged 
G55 equations at y = and y = 1/2 when T 55 = 0. 

This discussion then demonstrates the origin of the constraint: it is a consequence of 
requiring b = constant without a radion potential. With these assumptions the system is 
over-determined, and a fine tuning of the energy densities is required to maintain a static 
solution in the bulk. 
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2.2 The Effective Four Dimensional Action 



Below we derive the effective action for the four dimensional theory. We will demonstrate 
that the above explanation for the origin of the constraint equation can also be obtained from 
this 4D action. This should be expected, since without a stabilizing potential the radion 
modulus is massless and appears in the 4D effective theory. Consequently, the physics of 
maintaining a static modulus should also be obtained in the effective theory. 

To obtain the effective action, in the following we perturb the metric about the RS 
solutions as in Eqs. ( [2.10| ). We compute to 0(5 f) and to all order in b(t). Terms such as, 



e.g. dSa/dt ~ p ~ p 3//2 , are higher order in p, p* and are neglected. The effective Lagrangian 
is given by 

C eff = -M 3 f 1/2 dy^(R + A/M 3 ) 

J — 1/2 

-a 3 V r (b) + a 3 C Pl (a, ¥ Pl) , ...) + a 3 QiC TeV (Q b a, ^ eV \ ...) , (2.23) 

where V r {b) is the radion potential, and C denote the Lagrangians of the matter fields on 
the Planck and TeV branes, expressed in terms of the bare (not rescaled) fields, bare masses 
and induced metric. With our approximations the 4D effective action is then computed to 
be 

S eff = -^/^a 3 f(l-^^ + -o^^-^^^-/^a 3 K(&) 
k 777.Q J \ a z ab 4 tr J J 

+ J dta 3 C Pl + J dta 3 QlC TeV . (2.24) 



We have expressed the induced metric on the two branes as g^ v = Q% diag(l, 
Since the curvature scalar 



-a 2 , —a 2 , —a 2 ) 



~ a /a\ 2 . 

n (4 )(a =-6- -6 - 2.25 
a \aj 

integrating the second term by parts results in a more conventional looking action 

+S% + S% v (2.26) 

where S M are the matter actions on the two branes. 

It is now straightforward to compute the b equation of motion. The point is that due to 
the dependence of Qb on b, the presence of the matter on the branes generates a potential 
for b. To see this, we compute the variation of the above action with respect to b, noting 
that (we assume) S M depends on b only through the warp factor Thus the contribution 
of the matter fields to the b equation is 

SS" 6S M 6gV fz- ,„„S2i 



a is- a >^r-%™!' (2 - 27) 
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where = diag(l, —a 2 , —a 2 , —a 2 ), and T is the trace of the stress tensor[j] in terms of 
the bare fields and bare masses, and is equal to p — 3p for a perfect fluid. Thus the matter 
fields generate an effective potential for b that is 

Veff(b) = \(P*~ 3p* + (p- 3p)n$) , (2.28) 

where we have added for later convenience a b independent contribution from the Planck 
brane (which does not contribute to the b equation). 

For general p and p*, the minimum of this potential is at b —>■ oo. In fact with 6 = and 
V r = 0, the b equation is 



3S!j(a ioW V-^3 p -^4a , (2.29) 

2 k z \aj 




which simplifies to 




(3p-p)n* = ^ [-+(-) j . (2.30) 



We note that this is just the average of the G55 equation at the location of the TeV brane. 
Using Eqs. ( |2.14| ) and (|2.15| ) the above equation simplifies to 



(3p - p)n 4 b = - (sp* - p* + (3 P - P )ni) • (2.31) 

For generic p and p* the solution to this equation is fif, — > 0. Since Q b = e ~ m ° b ^ 2 , b — > 00. 
That is, the branes want to blow apart. There is another solution, however, and that is 
to allow a fine tune between and p. In fact, inspecting the above equation implies that 
3p* — p* = — fll (3p — p) is also a static solution. Combining this with the conservation of 
energy equations for p, p*, then implies the earlier constraint p* = — pVt 2 . 

So this demonstrates that the constraint p* = — pQl is quite general, as it does not 
depend on the details of the solutions in the bulk. More importantly, this effective action 
computation explicitly shows that the constraint directly follows from requiring that the 
radion modulus is static even when there is no stabilizing potential. From this perspective 
it is clear that with a radion potential p* and p will not be required to be correlated. For 
un correlated p and p*, the branes want to go off to infinity; this however, will be balanced 
by the restoring force from the radion potential. What was once a constraint equation for 
p and p*, in the presence of the radion potential becomes an equation determining the new 
equilibrium point: A6~0(p, p*). 



*Thus we see that b couples to the trace of the stress tensor, and has "dilaton" -like couplings. 
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3 The Solution to the Einstein Equations in the pres- 
ence of a Stabilizing Potential 

We have seen in the previous section that the constraint on the matter on the two branes 
is a consequence of the fact that without radius stabilization one has an extra light field 
in the spectrum, whose equations of motion are static only for correlated values of the 
matter densities. In this section, we will show that the conclusions reached above are in 
agreement with the detailed solutions of the 5 dimensional Einstein equations, once the 
radius is stabilized. To study this problem, we assume that a five dimensional potential U (b) 
has been generated in the 5D theory by some mechanism (for example see ||19||). Then the 
equations of motion in the bulk are given by 

G 00 = K 2 n 2 (A + 17(6)) , G u = -k 2 o 2 (A + U(b)) , 

G 5 = 0, G 55 = -K 2 b 2 (A + U(b)+bU'(b)), (3.1) 



with Gab given in |2.2| . To simplify the solution of the equations, we assume that the 
mass of the radion is very heavy, and that near it's minimum U is given approximately by 
U(b) = Mj (^ a ) 2 , where b is the stabilized value of the radius, and the parameter M\ is 
proportional to the radion mass m Ta di on - In order to understand what is happening to these 
equations we assume for a moment that Mj is the largest mass scale in the theory. Then the 
solution to the last equation in Eq. Q3.1J ) is simply given by b = bo, with no other constraint 
on a and n. With this solution we also find that U(b )=0. Thus in the presence of a heavy 
radion field the relevant Einstein equations are the 00 and %% components, with the radius 
fixed to be at the stable value bo- This already shows how the constraint is eliminated. One 
of the equations of motion, which played a vital role in establishing the correlation between 
the matter on the two branes is simply not appearing, it is automatically satisfied in the 
presence of the stable radius. 

The physics of this is that if the radion is very heavy, it can always adjust in an in- 
finitesimal way to satisfy this equation. As an analogy one can think of two charged spheres 
which gravitationally attract but electrically repel each other. There will be a static solution 
only, if the force from the electric charges exactly cancels the force from the gravitational 
attraction - thus there is a relation between the masses and charges of the two spheres. If we 
connect a spring with an extremely large spring constant between the two spheres, however, 
the only condition we will have is that the spring be at its own equilibrium point, and the 
other forces are canceled by infinitesimal changes in the length of the spring. 

Thus what we are left to show is that the first two equations, with the radion set to bo have 
solutions for arbitrary values of matter perturbations. To simplify the calculation, we perturb 
around the RS solution with cosmological constants 5V and 5V* instead of matter densities. 
The appropriate Hubble law for matter densities can be obtained by the substitution 5V — > p, 
5V* — > p* and H — > d (t)/a (t) (to lowest order in the perturbations). We solve the equation 
using the ansatz 

a(t,y) = e at e-M rn ° h °(l + 6a(y)) 1 n(t, y) = e -«'(l + 5a(y)), b = b . (3.2) 
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This way the 00 and the ii Einstein equations reduce to a single equation of the form 



HX = a{y)a"{y) + a\yf - 2m 2 b 2 a(y) 2 , (3.3) 

where a(y) = e~^ m ° bo (1 + 5a(y)). Plugging back this ansatz and expanding to first order in 
Sa(y) we obtain the equation 

8a" {y) - 4b Q m 8a'{y) = b 2 H 2 e 2 ^ mobo . (3.4) 

This equation is a linear inhomogeneous ordinary differential equation, which can be solved 
by the standard rules. The general solution is given by 

5a(y) = ^—( e *\y\ b «™« - i) _ _^( e 2|#om _ i) (3.5) 
4m 6 4mg 

where the overall constant has been fixed such that 5a(0) = 0. The remaining two constants, 
a and H 2 have to be fixed such that the jumps of this function at the two branes reproduce 
the matter perturbation that we are including. The result is given by 

and the value of the other constant a is given by 

a= f^ (6V.S% + WS%). (3.7) 

We note that Eq. ( |3.6| ) is the standard Hubble law formula with correct normalization for 
the physically observed energy density 8V* + Qq8V. The consistency of this solution requires 
8a -C 1, which by inspection implies that 8V* is much less than the intermediate scale 
figMp;, and that 8VQq <C (TeV) 4 . Thus the 4D effective theory should be valid when these 
conditions are satisfied, in particular if the temperature on the TeV brane is below the TeV 
scale. 

The lesson is as expected: radion stabilization removes the constraint between the matter 
fields on the two branes, there is a solution to Einstein's equations for any value of the 
perturbations, and the Hubble constant is given by the expression expected from the effective 
4D description of the theory. To show that our conclusions are generic for brane models, we 
derive a solution similar to the one presented in this section in Appendix A for the case of 
vanishing cosmological constants (the case considered in [E3]) in the presence of a stabilized 
radius. 



4 Radion Kinematics and Dynamics 

In this Section we concentrate on the radion field and point out some implications of the 4D 
effective action, Eq. (|2.24|) . First we discuss the cosmological consequences of ( |2.24j ), and 



then we examine the couplings of the radion to the SM fields and the consequences of these 
couplings. 
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4.1 Radion Cosmology 

Recall that for the metric (dropping the 5a, 5n and 5b perturbations of Section §, since they 
only contribute at 0((5f) 2 ) ) 

a(t,y) = a{t)e- mamy \ 

n(t,y) = e - m ° myl 

b(t,y) = b(t) (4.1) 

the 4D effective action is 

+S» + S% v . (4.2) 

From this equation we learned in Section 2.2 that without a stabilizing potential a constraint 
must be imposed on the different energy densities to maintain a constant radion modulus. 
In this Section we include a stabilizing potential V r . 

The presence of the interaction ab implies that in this coordinate basis the physical radion 
is mixed with the massless graviton. To separate the fields, it is convenient to perform a 
conformal transformation on the metric: 

a{t) = f{b{t))a(t) 
dt = /(&(*)) dt . (4.3) 

We want the action in the new basis to contain no cross terms ab. This fixes / to be 

/w = (bw) • (4 ' 4) 

where recall that Qq is the present-day value of Qf,. Note that for small perturbations of Q b 
away from its present-day value, f{b) = 1 — Q 2 ] m 5b/2 H — •. So for small perturbations, i.e. 
m 5b 1, the difference between the Einstein frame and the original frame is incredibly 
small, O(10~ 30 ). The difference between the two frames is then only important for large 
departures in f2& from its current value. 

Then the 4D effective action in this frame is : 



+ 



dia 3 f(b) 4 C Pl (a, ¥ Fl \ ...) + J dia 3 f(b) 4 niC TeV (n b a, ¥™\ ...). (4.5) 



In this action the radion appears as a scalar field with a non-trivial (6-dependent) kinetic 
term. In the absence of sources the radion potential is 

V r (b) = f{bfV r {b). (4.6) 
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In the presence of sources one can show similarly to the arguments in Section 2.2 that this 
is modified to 

V r , eff (b) = V r + ^ (p* - 3p* + (p - 3p)fi£) . (4.7) 

It is also convenient to define a new scale 

A w = Q M Pl . (4.8) 

Then Ky/ ~ 0(TeV). Inspecting the above Einstein frame Lagrangian, we see that for small 
perturbations away from f2& = Q , the canonically normalized radion cf) is 





™ob(t) = k -^(1 - ng) ~ a . (4.9) 



In the following the corrections are neglected, so the second expression on the RHS is 
used. 

In this frame the Hubble law has a simple expression and interpretation: (This can also 
be obtained from performing the conformal rescaling in Eq. (|2.14 )) 




a 8ttGj\[ ( r . nX A, s 13,, 

— = —5— /(&) p* + pvis) + -;^-7r- ( m ob) 
a 1 3 \ 4 o7rG tv 

= ^ (^0 2 + VM) + f(b)\ P * + Pms )) ■ (4.10) 

Here and throughout Gn is the present-day value of Newton's constant, and the second 
expression is only valid for b close to bo. We see from this equation that as long as the 
energy density in the radion does not dominate the energy density in p vis = pQf, in the 
Einstein frame the universe expands as in the usual cosmology. Since by assumption the 
oscillations are small, in the original frame we also have FRW expansion since then a(t) = 
a(t) x (1 + small osc). 

In the absence of sources (p = p* = 0) we require that the potential truly stabilizes the 
radion. That is, we require b = is a solution. Inspecting the b equation we see that a static 
solution in this case is possible only if V =0. If we also require that a is static in this limit 
then the a equation of motion implies that V = 0. Since / 7^ and /' 7^ 0, at the local 
minimum we must have V = and V = 0. With this knowledge, the radion mass obtained 
from the above action is then given by 

2 _ 2 W(&o) q n2\2 (AU) 



Near the local minimum it is then convenient to expand V r as 

Vr(b) = ^f^V^M 




(4.12) 
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For the Goldberger-Wise stabilizing mechanism, V" ~ Mj> z x O(TeV) 4 (see Appendix B), so 
from the above formula m r ~ O(TeV). This is not surprising, since a bulk scalar field with 
O(Mpi) bare mass appears in the 4D effective theory as a Kaluza-Klein tower of scalars, 
where the lightest field has an O(TeV) mass. 

In the following the action for the canonically normalized radion <ft is then truncated at 

Sradion = \j ^Xy/g ((<90) 2 - 77^(0 - O ) 2 ) • (4.13) 

This will be the appropriate action for describing small displacements of the radion from its 
minimum. 

In Section ^| it was argued that the special constraint between the energy densities on 
the two branes disappears once the radion is stabilized. Using the Einstein frame action, 
Eq.( f4.5|) , the b equation is 

3 2, f'b Jb m bb 2 l + Q 2 \ -, , A1A . 



16nG N (l-n 2 b ) 2 \b ab 2 b 2 l-n 2 )~ r ' eff 

This equation can be used to recover the earlier constraint p* = — Vt 2 p when V = and 
b = is imposed. Once this constraint is not satisfied, it is clear from inspecting the 
above equation that the restoring force term V will balance the expansion force term. More 
concretely, plugging V r , Eq. (|4.12|) , into the above equation and neglecting the b and b terms 
gives 

Ab _ 1 1 - ttl (p vis - 3p vis ) 

b ~ 3 m b m 2 r ttlM 2 Pl 1 ' ' 

where Ab is the distance between the minimum of the effective potential with and without 
matter, and we assume that the radion was at its minimum at the start of matter-domination 
(MD). Also p vis = f2 4 p, etc. is the physically measured energy density on the TeV brane. To 
obtain this formula it was also assumed that the Planck brane is empty (p* = 0). To include 
p* 7^ 0, in the above formula replace p vis — > p vis + f2 2 p*, etc. The consistency of neglecting 
the b and b terms requires Ab/b <^ 1, for then b ~ p ~ p 3//2 is suppressed compared to the 
Ab term. 

From Eq. ( p^D it is seen that there is no shift in b during a radiation dominated era 



(w = 1/3) from the dominant component.^ Hence the physical implications of this shift 
are much smaller than what one would have naively thought. During this era, however, the 
energy density contained a small component of non-relativistic matter, which will cause b 
to shift. Since both Newton's constant and SM particle masses depend on the vev of the 
radion, a substantial shift in these quantities during any era after BBN would lead to a BBN 
somewhat different from the standard BBN cosmology. The success of standard BBN then 
implies that these changes must be small, and this is found to be the case here. To estimate 



*This is because the stress tensor of radiation is traceless. In what follows we neglect the quantum 
corrections to Tr T. 
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the shift since the start of BBN, we note that the energy density in non-relativistic matter 
just before the start of BBN was p^n ~ (T BBN /T ) 3 p Cj0 ~ 10 20 eV 4 , where we have used 
Tbbn ~ 10 MeV, p Ci o ~ 1CT 10 eV 4 . Substituting p vis — 3p V i S ~ pnr into Eq. (|4.15|) gives 



6 °N . o2_ jfi. m-30 w 1n -28 ^TeV\ 2 

Gat 



nimoSb ~ 1(T JU x 1(T 28 — . (4.16) 

V m r J 



We also require that the weak scale has not shifted by more than (9(10%) since BBN, which 
implies Ab/b < (9(l/m &o) • From Eq. ( (4.15[ ) this implies 



ao(eV)' 



m r > 10 iU ~ O(10- 2 ) eV. (4.17) 

lev 

Of course, it is obvious from the discussion in Section 4.2 that accelerator experiments 
provide a much stronger constraint on the radion mass than the above result. 

One general drawback of the RS model seems to be that there are two fine-tunings re- 
quired to obtain a static solution (see Eq. (|2.6|) ). One of these tunings is clearly equivalent 
to the usual cosmological constant problem, however, it seems that there is a second equally 
bad tuning appearing, making the model less attractive. We will argue that (as already 
emphasized by |19, (43| ) the second tuning was only a consequence of looking for a static 



solution without a stabilized radion, and disappears in the presence of a stabilization mecha- 
nism. To see this, assume that the radius is stabilized, and we perturb the brane tensions as 
V — ► and V* — > V*+5V* so that the fine-tune relations Eq. fl2.6| ) no longer hold. From 

the previous discussions about stabilizing the radion, we find that the first perturbation gen- 
erates an effective potential for b: SV e ff ~ SVCl^. This perturbation has two effects: to shift 
b, and to increase the 4D cosmological constant. The 4D cosmological constant is canceled 
by appropriately choosing 5V*, so one fine-tuning equivalent to the 4D cosmological constant 
problem remains. But if 5V is small compared to both V and the typical mass scale in the 
stabilizing potential, then from Eq. ( 4.15|) we find that the shift in b is tiny. Thus once the 



radion is stabilized, only one fine-tuning of O(Mpf) is required (equivalent to canceling the 
4D cosmological constant). The other tuning required in order to get the correct hierarchy 
of scales is clearly only of (9(10 _1 — 10~ 2 ). The conclusion is that there is only one severe 
fine-tuning in the RS model, equivalent to the cosmological constant problem. 

We next argue that there is no moduli problem associated with the radion. A potential 
moduli problem can occur from two sources. The first is that coherent oscillations of the 
radion can overdose the universe well-before BBN, ruining the success of BBN. In the next 
section we will see that after EWSB, the radion has renormalizable interactions with the SM 
fields. So after EWSB it promptly decays, long before the start of BBN. So by the start of 
BBN there is no energy in the coherent oscillations of the radion. The second concern is that 
the above shift in b may introduce a correction to the Hubble Law which would modify the 
time evolution of our universe. This too could spoil the success of standard BBN cosmology. 
We shall find that the correction is sufficiently small that this is not a concern. 

The shift in b will add energy to the radion potential. As the universe cools this energy 
will not red shift like matter, since it is not a coherent oscillation. Rather, as the universe 
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cools b will follow its instantaneous minimum since the cooling is adiabatic (for m r 3> H). 
The energy stored in the radion potential then appears as a p 2 correction to the Hubble law. 
This could affect the time evolution of the universe. For a general equation of state, we find 
from the above expressions for Ab and V r that 

V r (Ab) (3p vis - p vis ) 2 p 2 NR 



Pvis m r-h-wPvis ^r^wPvis 

Using the value of p NR ~ 10 20 eV 4 when T ~ 10 MeV, then the LHS is < 1 if m r > 10" 6 eV 
which is weaker than the collider constraints. So there is no moduli problem from the shift 
in the radion. 

We conclude this section by arguing that a stabilizing potential is in fact required in 
order that the SM particle masses have not changed significantly since the start of MD. (Of 
course direct accelerator searches require that the radion is heavier than (9(10-100) GeV.) 
The point is that without a stabilizing potential, the source terms on the TeV wall generate 
a potential for b, which wants the walls to separate. We can estimate the rate of expansion 
for b by setting the stabilizing potential to zero in ( f4.14p (which does not mean that the full 
V e ff vanishes). Inspecting that equation one can find that the ab term never dominates the 
expansion, thus the radion is not slowly rolling. One finds that the b and b terms in ( |4.14 ) 
are comparable, and this results in an expansion rate 

I ~ (4.19) 
b A w 

for b close to bo- Thus the fractional change in the radius will be significant (O(l)) over time 
scales much shorter than the Hubble time, which results in significant changes in the SM 
particle masses during MD, completely changing the predictions for BBN. Thus the success 
of BBN cosmology requires that the radion is stabilized before the onset of BBN. 



4.2 Couplings of the Radion to SM Particles 

The radion field interacts with SM particles, and the size of these interactions determines 
whether the radion can decay quickly enough to avoid overclosure by the start of BBN. We 
shall see that the radion has renormalizable and and 1/TeV suppressed interactions with the 
SM fields. 

In the RS scenario, the radion has renormalizable interactions with SM fields after elec- 
troweak symmetry breaking (EWSB). This is because the mass scales on the TeV brane 
depend on the radion modulus by 

<h>=v = v e . (4.20) 

But since this warping contains the radion modulus, this dependence introduces renormal- 
izable couplings of the radion to matter once EWSB occurs. For example, in the Yukawa 
interaction we expand the Higgs vev as in Eq. (|4.20|) to obtain 

- - 1 v — 

Xijh^j -> XijVi/;^ - X i:j /gA^ ^M'j + • • • • (4-21) 
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So if present, this results in renormalizable interactions between the radion and the SM 
particles, scaled only by v/A w in addition to the usual factors of Yukawa couplings or gauge 
couplings. Since the radion mixes with both the neutral Higgs and with the massless graviton 
(though the latter is negligible compared to the radion-Higgs mixing), we must work in a 
basis where the radion does not mix with these fields. 

Next the action for the radion and Higgs interactions is derived, and the physical eigen- 
states and their couplings to SM fields are identified. We parameterize the metric on the 
TeV brane as 



and in the bulk as 



g^u(x,y = 1/2) = tt 2 b (x)g^ v (x) , 

n b {x) = e - mob{x)/2 

fafay) = e- 2m ° b ^g» u (x) , 
gssfav) = ~b 2 (x) . 



(4.22) 
(4.23) 
(4.24) 

(4.25) 



This expression then includes the fluctuations in the radion and the zero-mode graviton. 
Inserting this into the 5D action and integrating over the extra dimension results in the 
following 4D effective action 

Seff = J d 4 x^ (-l^^-n i4) (g) + ^n 2 b (db) 2 - Vribfj + S TeV + S Pl .(4.26) 

From this action we find that for small fluctuations about b = b , the canonically normalized 
radion is 

mMx) = fJM (4 , 7) 

which is the same as Eq. ( f4.9| ) for small Qq. (The small difference between this and Eq. ( |4.9| ) 
is due to the latter normalization being defined in the Einstein frame.) The SM action on 
the TeV brane is 

S Te v = J d A xJ~g [g^D^D v E - V(H) + i^e^D^ - \ ij E^ j + h.c.) 

= J d A x^n A b {nfg^D^DvH - V(H) + in^e^D^ - X^H^ + h.c.) , 

(4.28) 

where for simplicity only the fermion and Higgs interactions are included. It is straightfor- 
ward to include other interactions, for example, the Yang-Mills kinetic terms. We assume a 
Higgs potential of the form 

V(H) = X(H^H - \vl) 2 . (4.29) 
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To obtain canonical normalization of the kinetic terms in Eg. Q4.28 ) we rescale H — > Q b 1 H 
and ip — > Vt^^il) . The effect of this is to renormalize Vq — ► fi&t>, and to introduce higher 
dimension operators between the radion and the SM fields. That is, the new SM action is 

S ef f = J d 4 x^g (g^D^D v H - V(H, 0) + i&e^D^ - X^H^ + h.c. 
+ S radion + / d^ ^ DH + h . c) + (J| 




+ J d 4 x^g ^—^—ty^d^ j+S Pl . (4.30) 

Note that in the kinetic terms the radion has dimension 5 and dimension 6 interactions with 
fermions and bosons. The Higgs potential after the rescaling is 



V{H,<j>) = A (h*H- ^o e_V ^^) • ( 4 - 31 ) 

Note that this confirms the intuition that the source of the radion coupling to SM particles is 
the conformal breaking sector. This potential generically results in mass mixing between the 
Higgs and the radion. Inspecting the final action, Eq. ( |4.30[ ), we also see that after EWSB 
there is kinetic mixing between the radion and Higgs. So in order to verify the claim that 
after EWSB the radion has renormalizable interactions with the SM fields, we must identify 
the coupling of the SM fields to the physical eigenstates made out of H° and 0, and in 
particular verify that the gauge eigenstate H° is not a mass eigenstate. 
Expanding H° = [v + h + iHi) /y/2, the and H° equations are 



and 



= m 2 r {<j) - O ) + X-fv e ^ A w ( v 2 -v\e v 3 A w ) (4.32) 



= \v(v 2 - v\e v? 4v ) , (4.33) 



where the scaling factor 7 is defined by 

Note that 7 is fixed by the precise value of Qq- So the vacuum that breaks EWS is 

< > = 00 

v 2 = t) 2 e V3 % . (4.35) 

This is physically reasonable, since the radion and Higgs potential is the sum of two positive 
definite quantities, which has a local minimum when each term separately vanishes. This 
also demonstrates that at tree-level the fermion and gauge boson mass spectrum is the same 
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here as in the SM, since the vev of <p only determines v, and otherwise does not contribute 
to any of those masses. 

In the original gauge basis the mass matrix for (0, h) is 



/ ml + 2 7 2 Av 2 2 7 Aw 2 \ 
^ 2 7 Aw 2 2\v 2 ) ■ 

In this basis, however, the kinetic terms are not diagonal, but are given by 



(4.36) 



Uin = \{l +l 2 )(d<P) 2 + ±(dh) 2 + 1 d<pdh. (4.37) 

The kinetic mixing is undone by the field redefinition^ (f> — </>o = and h = h — j<f>. In the 
(0, h) basis the mass matrix is also diagonalized. So (<f), h) are also mass eigenstates, with 
mass squared m 2 and 2Av 2 , respectively [[ Note that these are unchanged from their naive 
expectation. The relation between the gauge and physical basis is then 

<\> = 0o + 4> (4-38) 

and 

h = h-~/4> . (4.39) 

Thus the couplings of the radion to SM fields are similar to the neutral Higgs, and may 
be important for collider phenomenology since 7 may not be small. We point out that the 
linear couplings of the radion to SM fields may also be obtained from : 

= = ^ = 7-T • (4.40) 

So to 0(4>) the radion couples to the trace of the stress-tensor. To this order, this is equivalent 
to the coupling of the radion to the SM in the large extra dimension scenario, except that 
here the suppression is TeV rather than Mpi. 

This discussion then demonstrates that at the EW scale the theory contains two neutral, 
CP even, spin-0 particles coupled to the SM fermions and gauge bosons. The tree-level 
mass parameters are the same, however, as in the SM, since the vev of the radion does not 
break electroweak symmetry. After EWSB the radion has both mass and kinetic mixing 
with the neutral Higgs, which after diagonalization results in renormalizable interactions 
between the physical radion and the SM fields. These interactions are scaled by vj \^/6A w ), 
so the radion is more weakly (or more strongly!) coupled to SM fields than the neutral 
Higgs. We also point out that the radion also couples to the kinetic energy terms, and these 

' Note, that the radion also has kinetic mixing with the graviton. Due to the weak dependence of Mpi on 
b, however, this mixing is proportional to £7^, and is thus negligible compared to the Higgs mixing considered 
here. 



J This is also seen if in Eq. ( 4.28 ) we rescale the fields by the vev of fib instead of the full quantum fib. 
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operators become important at EW energies. It is then interesting to determine the current 
experimental limits on the radion mass, and to what extent future collider experiments can 
distinguish between the radion and neutral Higgs boson. We note that these important 
experimental issues are completely determined by the precise value of Qq. 

The dominant decay of the radion is to it, ZZ or WW if kinematically allowed, otherwise 
to bb. At e + e~ colliders it can be produced in association with Z, so its collider signatures 
are similar to the neutral Higgs, but with suppression or enhancement of the rate. The Higgs 
may also decay to two radions (if kinematically allowed), and may be competitive with its 
SM decay modes. 

The cosmological implications of this is that the radion can decay to SM particles before 
the start of BBN. Of course, a requirement on early cosmology for this to occur is that by 
the start of BBN the radion was close to its present-day value. Then the lifetime of the 
radion was small enough, and the energy density in coherent oscillations of the radion was 
be transfered to radiation by the start of BBN. 

We also note that if matter is present on the Planck brane, the coupling of the radion to 
that brane is 0(1/Mpi) or smaller (see Appendix E). Consequently, the radion has a very 
small branching fraction to decay on the Planck brane. This means that in the early universe 
the coherent oscillations of the radion would predominantly dump their energy into the TeV 
brane, and only transfer a small fraction 0((TeV/Mp;)) 2 to the Planck brane. It also implies 
that during the time that b was close to its present-day value, matter on the Planck brane 
was not in thermal equilibrium with the matter on the TeV brane. 

We conclude this section by noting, that (just as for the case of the KK gravitons) the 
radion of the RS model has very different properties from the radion of the large extra 
dimensions of ||]. In the case of large extra dimensions, the lower bound on the mass of 
the radion is 10~ 3 eV from the gravitational force measurements. In these models there is 
also an upper bound on the radion mass, which is obtained by requiring that the curvature 
radius is smaller than the physical size of the extra dimensions. This upper bound varies 
between 10 -2 eV for two extra dimensions to 20 MeV for six extra dimensions. Thus in these 
models the radion is typically very light, with mass between 10~ 3 eV and 20 MeV, while its 
couplings to the SM fields are suppressed by Mpi. In the RS model however, the radion has 
unsuppressed couplings to the SM fields, which result in a lower bound of order 10 — 100 GeV 
for the radion mass, while there is no upper bound from the curvature constraint. Thus one 
can see, that the radion in these two models has very different properties, and they should 
be easy to distinguish experimentally. 

5 Dark Matter on the Planck Brane or in the Bulk? 

The expansion rate of the universe is driven by three sources: matter living on the TeV 
brane, matter living on the Planck brane, and matter living in the bulk: 




(5.1) 
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where p bu i k = b J\Q{y) 4 p bu i k {y)dy, with p bu i k {y) being the five- dimensional bulk energy 

2 

density. (Recall that for late-cosmology it is a good approximation to set V r = and 6 = 
in the Hubble law. Also the time dependence in b due to the shift in b contributes 0(p 2 ) to 
H 2 which is negligible for late-cosmology). 

We see from Eq. (|5.1|) that the appearance of and p bu[k in the expansion rate is 
potentially dangerous. This is because the natural scale on the Planck brane is 0(M Pl ), 
and so this energy density could easily overdose the universe. But we will argue below 
that natural mechanisms exist to suppress the energy density on the Planck wall well-below 
its natural value. But of course a careful computation within a specific model for early 
cosmology is required to establish that these suppressions are indeed suffucient. 

Since the natural scale on the Planck brane is M pl , we see from the Hubble Law that 
there are two different energy scales appearing in the effective 4D theory, which is at first 
surprising. This suggests that observers living on the TeV brane will measure Planck-mass 
particles living on the Planck brane to really have Planckian masses, even though the cutoff 
on the TeV brane is O(TeV). In fact, this is confirmed in Appendix D by a computation of 
the Newtonian force law between a particle on the Planck brane and a particle on the TeV 
brane. We find that the Newtonian force law between particles of bare mass ni\ and rri2 
located at y\ and y 2 along the extra dimension is given by 

F N = G N m 1 Q(yi)m 2 Q(y2)/r . (5.2) 

So the Newtonian force law and Hubble Law for matter living on different branes are con- 
sistent. 

Thus one can see from (|5.1|) , that there are two new candidates for dark matter: matter 
in the bulk and matter on the Planck brane. The radion is not a candidate for dark matter, 
since as we have seen in the previous section, it has a short lifetime (of the order M~^ ak ), 
and is therefore not abundant today. 

Matter in the bulk can naturally be dark matter, since its natural mass scale is 0(TeV)[], 
and it has interactions with the SM fields suppressed by the TeV scale. Thus it could have 
been in thermal equilibrium with the TeV brane when the temperature of the Universe was 
0(TeV). If these bulk fields are unstable, then they will predominantly decay to SM fields, 
rather than matter living on the Planck brane, since there the couplings are suppressed by 
Mpi. However, if those bulk fields are stable, and since it is natural for their couplings to 
SM particles to be TeV suppressed (see Appendix E), their annihilation cross section may 
turn out to lead to interesting relic densities today, which may serve as dark matter. 

Another interesting possibility is to have matter on the Planck brane serve as dark mat- 
ter. The natural mass scale for particles on the Planck brane is either Mpi or (massless). 
In either case one might worry that matter on the Planck brane would overdose the Uni- 
verse. However, matter there is very weakly coupled to bulk fields, since its couplings are 
suppressed due to the small wavefunction of the bulk fields at the Planck brane, in addition 



*A 5D scalar with O(Mpi) bare mass appears in the 4D effective theory as a Kaluza-Klein tower of scalars. 
The lightest mode has a mass of 0(TeV)@. 
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to the suppression of Mpi for non-renormalizable operators. We illustrate this in Appendix 
E, where we calculate the couplings of bulk scalars to the Planck (and the TeV) branes. 
Therefore matter on the Planck brane has presumably never been in thermal equilibrium 
with our brane after inflation. Even if massive matter there was in a thermal equilibrium 
with the TeV brane (or if the temperatures on the two branes were once comparable), its 
density is suppressed by a Boltzmann factor, and is therefore greatly suppressed. Such a 
suppression would however not happen for massless fields on the Planck brane (unless the 
temperature itself on the Planck brane is for some reason much smaller than on the TeV 
brane), thus the presence of such massless fields is strongly disfavored, which is a constraint 
for model building. It would be worthwhile to investigate whether non-thermal production 
of Planck brane matter could result in interesting relic densities [42 and thus serve as dark 
matter. 



6 Conclusions 

In this paper we have considered the effect of a stabilization mechanism for the radion on 
the cosmology of the RS model. We have found that the previously discovered unconven- 
tional cosmologies result from the fact that the radion was not stabilized in the basic RS 
scenario. The constraint between the matter on the two branes was a consequence of trying 
to find a static solution to the radion equations of motion without actually stabilizing it. 
Once a stabilization mechanism is added, the constraint disappears, and the ordinary four 
dimensional FRW Universe is recovered at low temperatures if the radion is stabilized with 
a mass of the order of its natural scale M wea ^. 

In this paper we have also pointed out that after electroweak symmetry breaking the 
radion has renormalizable couplings to matter on the TeV brane. These interactions are 
uniquely determined by the precise value of the warp factor. The production and detection 
of the radion are then similar to the detection of the SM Higgs, but with enhanced or 
suppressed rates. An important consequence of these interactions for cosmology is that 
there is no radion moduli problem, as the radion decays into SM particles long before the 
start of BBN. We have also speculated that matter in the bulk or on the Planck brane are 
new candidates for dark matter. 



Note Added 



While this paper was being concluded, several papers related to this subject appeared |44, 
45 , 46] , however none of these papers consider the central issue of our paper, the effect 
of radion stabilization on the cosmology of the negative tension brane. Ref. derives 
the four dimensional gravity equations in a covariant formalism, and concludes similarly to 
55] , |56| : the expansion of the positive tension brane is conventional, however, there is a sign 

Refs. 



difference for the case of the negative tension brane. 



45], consider the cosmology 



of the positive tension brane. Ref. E5J concludes that after inflation one can always end up 
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with the "correct" expanding solution found in |35[. Ref. [[RJ] finds the exact solutions in 
the bulk to the Einstein equations with matter on the branes and background cosmological 
constants. The solution presented in Appendix A of this paper is a perturbative expansion 
of the exact solution of |46|, obtained independently from f46|. 
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Appendix 

A The Solution without a Stabilized Radius 

In this Appendix we find an approximate solution to Einstein's equations in the RS model 
with matter on the branes for the case when the radius is not stabilized. The main purpose 
of this is to illustrate the generic feature of all such solutions of brane models without a 
stable radius: there appears to be a constraint between matter on the two branes. This is 
similar to the relations found in [33[ 34, 38| . It has been suggested |33f that such relations 



are of topological origin and are very generic to brane world scenarios. In the Sections [| and 
|3] we argue that the real origin of this constraint is the requirement of a stable radius even 
in the absence of a stabilizing potential. Otherwise the extra dimension itself will expand or 
contract. The presence of the stabilizing potential then prevents the expansion of the extra 
dimension, and a phenomenologically acceptable cosmology follows. 
In order to find a solution we assume the form of the metric to be 

ds 2 = n(y,t) 2 dt 2 - a(y,t) 2 (dxl + dx\ + dx\) -b 2 dy 2 , (A.l) 
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that is we assume we have a time-independent constant radion b even though there is no 
stabilization mechanism. We also assume that the matter density on the branes is small 
compared to the brane tensions p, p* <C V. Then one can find a solution as a perturbative 
series in p around the RS solution: 

oo 

a(y,t) = a (t)e-^ bom °(l + Y,pMfi(y)), 

1=1 

oo 

n(y,t) =e-^ m °(l + J2p*(ty gi (y)). (A.2) 



1=1 



One could start with a more general ansatz, where instead of explicitly expanding in p*(t) 
we expand in arbitrary functions of time. However, the "jump equations" discussed below 
will immediately tell us that this function of time is proportional to p*(t). 

Since p* <C V, we will keep only terms that are linear in p*, thus our ansatz for the 
solution will be 

a(y,t) = a (t)e~M b ° m °(l + p*(t)f(y)), 

n{y,t) = e'\y\ b ^{l + p^t)g{y)). (A.3) 

We know from J33|, that the jumps of a and n are related to the density and pressure in the 
following way: 

MM>]- 4(V + P ,(t)), 



a(0,t)b 3 

KM £(-V + p m (t)(2 + Zw.)), (A.4) 



n(0,t)&o 3 

where [h] denotes the jump in the function [h(0, t)] = (h(0 + e) — h(0 — e)). From these 
equations we learn what the jumps of / and g have to be: 

[/'(0)] = 

[ff'(O)] = k%(w, + (A.5) 

The G 05 equation tells us that f(y) and g(y) have to be proportional to each other (up to a 
constant): 

p* + 3%*( l + ty^M] =0, (A.6) 

which implies that g'{y)/ f'(y) = const., and using the values of the jumps of the derivatives 
together with the fact that we have an S l /Z 2 orbifold implies that the G 05 equation just 
reduces to the usual conservation of energy equation of the form 

p* + 3— (p* +p*) = . (A.7) 
a 
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Similarly, at the other brane we have 



p + 3^(p + p) = 0. (A.8) 

Plugging the ansatz into the other equations then implies that 

f(v) = --^(e 2 '^ b » - 1), 
12m 

= k 2 (2 + 3w*) 2lylmQbQ _ 
yvy; 12m 

The constant in / and g have been determined such, that in the limit ttlq — > the solution 
exactly reproduces the bulk solution presented in [33|. In order to satisfy Einstein's equations 
the function a Q (t) has to satisfy the following Friedmann-type equations: 

a \ 2 1 2 

— J = o K m 0P*> 

a J 6 

a \ 2 a K 2 m 

— H = — — (p* - 3p*J. (A.10) 

ao / ao o 

From (|A.9|) one can see that the Hubble constant at different points along the extra dimension 
is given by 



(^y) a = R2{y) = H2{0) i 1 + 2^ (e%|mo '° ~ ^ +p ^ + ' (A - n) 

which for p # Mp ; is a very slowly varying function of y. This can be seen by comparing 
H 2 (\) and H 2 (0): 

i K 2 m b 

H\) - H 2 (0) ~ ^ 2 (0)^^(P* + P.). (A.12) 

If e m ° fe o ~ 10 30 , and ^ ~ then H 2 (D~H 2 (0) ~ g 2 (0) (^^i ■ Thus for temperatures 

below 10 10 GeV on the Planck brane the two branes will expand together, but the expansion 
rate is completely fixed by the matter on the Planck brane. This shows, that for these 

temperatures there should be a sensible effective four dimensional theory describing the 

i 

evolution of the Universe, since a(0) ~ a(^) ~ (a) = f\ a(y)dy. Below we will show that 

?. 

the effective theory calculations do indeed reproduce (|A.10|) . 

From the point of view of the negative tension brane one can understand the previous 
Hubble Law from the following consideration. The above solution depends only on p* and 
p*, and not on p or p. This means then that the jump equation at y = 1/2 completely 
determines the energy density on the negative tension brane. One finds that the condition 
is simply given by 

p*(t) = -p(t)n 2 , P *(t) = -p(t)n 2 . (A.i3) 



2 
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This explains the results found in |36|j|*|: the expansion equation has the wrong sign if 
one assumes p > 0. In |35| p > was assumed in order to get a sensible phenomenology on 
the TeV brane. However, we see that the constraint (|A.13|) implies that for p > one needs 
to have p* < 0. But from (|A.10|) we see that p* cannot be negative. Therefore, in this case 



p on our brane must be negative. This is completely due to the constraint ( |A.13| ), and not 
a consequence of the breakdown of the effective 4D theory. 

One can check this result by explicitly calculating the effective 4D action for this setup, 
using the insight from the solution. The full 5D action is given by 

- / d " x V9(^ + A) + / d 4 x^f d C Pl + J d 4 xJgfjl£ TeV . (A.14) 

To calculate the effective 4D action one needs to integrate over the extra dimension. We do 
this substituting for the metric the ansatz 

a(t,y) = a (t)(l-^(e»"-l)), 
n(t,y) = (1 + ^g±^l (e %l-o 6 o _ (A . 15) 

however without assuming that ao and p satisfy the above Friedmann equations. The result 
of integrating over the fifth coordinate y to linear order in p and p* is 

S eff = [ dta (t) 3 3{1 ; Ql) ( V + *)+ [ dta (tfC Pl + [ dta (tfQ^TeV 

J K TTLq V \CLq/ CLq J J J 

= J dtalit) (-±M Pl K^ + Cpi + ng/w) • (A. 16) 

This is exactly what we expect: the action for the expanding universe with the correct Planck 
constant, and with total energy density ptotai = P* + pfig- As expected, the energy density 
on the "TeV" brane contributes pQ$ to the expansion rate. The expansion rate obtained 
from this effective action also agrees with the expansion rate for the induced metric in the 
five-dimensional theory. In fact, the effective theory gives 

tf 2 = ^p(p*+P^) , (A.17) 

which after substituting the constraint Eq. ( |A.13|) , and the relation between Mpi and the 
five-dimensional parameters, Eq. ( |2.7| ), identically agrees with the expansion rate Eq. ( |A.1U| ) 
in the five-dimensional theory. 

Thus we can see that the effective theory picture is as expected in agreement with the 
detailed form of the solutions to the Einstein equations: there is no breakdown of the effective 

*Note, that a factor of J1q has been erroneously omitted in (14) of |35| and in (5) of |36| in the case of p_. 
The reason for this extra factor is that from the point of view of the negative tension brane the fundamental 
scale of gravity is of the order of MpjQg ~ 1 TeV. 
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theory picture since the two branes expand together. The only puzzle that remains is why 
one needs to have a relation between the matter fields on the two branes of the sort in 
Eq. ( |A.13| ). As we will see in Sections 2 and 3, this is a consequence of the fact that we 
have not stabilized the radius. It simply tells us, that if we put in matter that does not 
obey this constraint the extra dimension will want to expand itself. However, if the extra 
dimension is stabilized, no such constraint should exist and the cosmological expansion is 
given by the conventional Friedmann equations in the 4D effective theory with no constraint 
on the magnitude of the matter one can add to the system. 



B The Case of Vanishing Cosmological Constants 

In order to show that our results are generic for any brane model, we will convince ourselves 
by using arguments similar to those in Section 4 that the ordinary Friedmann equations 
are recovered for the case of vanishing background cosmological constants as well (the case 
considered in f33|), once the radion is stabilized. To see this in detail, we again look for the 
solutions to the 00, ii and 05 components of the Einstein equations only, and just like in 
Sec. |H we do not require that that the 55 component is satisfied, since this will be solved 
by adjusting the radion, which we once again assume to be very heavy. We assume that we 
include matter density p(t) on "our brane", and matter density p*(t) on the other brane. 
(Note, that since the radius is stabilized, the presence of a second brane in the case of 
vanishing background cosmological constants is not even necessary. One can obtain all the 
relevant formulae for this case by just setting p* = 0.) We look for an approximate solution 
(i.e. valid to 0(p*,p)) on the S 1 /Z 2 orbifold in the form 

a(t,y) = a (t)(l + ap*{t){y - + (3p{t)y 2 ), 

n(t,y) = (1 + 7 p.(t)(y " \f + Ap(% 2 ), 

b{t,y) = bo{l + 6b), (B.l) 

for < y < |, and the solutions for the other regions are obtained by reflecting around y = 
or y = |. In this ansatz a, /?,7 and A are constants. From the jump equations we obtain 
that 

_ n 2 b (2 + 3w)K 2 b A (2 + 3w*)n 2 b 
a = (3= ^ 7 = 6 ' A = 6 ' (R2) 

where p = wp and p* = u^p*. One can easily check that with this choice of constants the 
relevant components of the Einstein equation are satisfied to first order in K 2 bop which is 
assumed to be small, and ao satisfies the following Friedmann equations: 

d \ 2 k 2 
a J 3b 

+ 2^ = ~(wp + w.p.). (B.3) 
a J a b 
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These are the correct 4D FRW equations with the correct normalization, since H 2 oc p, rather 
than H 2 oc p 2 , and as Mp t = bo/n 2 is the relation between the 4D and 5D parameters. The 
size of 5b is obtained from inspecting the G55 equation and is readily seen to be 0(p 2 , pi, pp*). 
Once again, we find that after stabilization of the radion the ordinary four dimensional 
Friedmann equations are recovered, with no constraint on what kind of matter one can 
include on the branes. Thus the ordinary FRW cosmology is recovered which results in 
ordinary BBN. 



C Radion Mass in the Goldberger-Wise Stabilization 
Mechanism 

Here we compute the radion mass using the Goldberger-Wise (GW) mechanism for radion 
stabilization [[HJ. It is found that the radion is naturally of O(TeV). 

In their mechanism for generating a stabilizing potential, a bulk scalar field $ is intro- 
duced which has a 5D action 

S 9 = ~ [ <£>x\[g~ (G AB d A ^d B ^ - m|$ 2 ) , (C.l) 



2 

where Gab is the full 5D metric which includes the RS warp factor. The bulk mass m| is 
assumed to be somewhat smaller than the Planck scale. So 

i = M (C.2) 

is a small quantity. The bulk scalar also contains two potentials on the Planck and TeV 
brane, 

S Pl = J d A x^\ h {& - v 2 f (C.3) 



and 



S TeV = J d'x^gX^ 2 - v 2 v ) 2 . (C.4) 



Note that v v and Vh have mass dimension 3/2. GW solve the equation of motion for $, and 
find that it has a non-trivial dependence in the bulk. In the large A„,h limit, GW find that 
$(0) = Vh and $(1/2) = Vh, so that in this limit the brane-potentials do not contribute 
to the 4D energy. However, $ contains both potential and (gradient) kinetic energy in the 
bulk, so that the energy density in the bulk is inhomogeneous. The radion potential is then 
obtained by substituting the classical solution for $ into the above action and integrating 
over the extra dimension. GW find that the resulting 4D effective potential for the radion 
is given by 



V r (b) = Am e- 2mob (v v - v h e- e ' m ° b/2 ) 2 (1 + -) - e'm v h e- {4+t ' )m ° b/2 (2v v - v h e- € ' m ° b/2 ) (.C.5) 



4 
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In their computation GW dropped terms of 0(e' 2 ) in V r , so the following computation of the 
radion mass is only accurate to this order. Then 

K-'W / . .. -e'mnb/2\ I n/_. -e'mn6/2\ , 6 ' 



- (y v - v h e~^ b ' 2 ) 2(v v - v h e~^ b ' 2 ) + -(v v - 3v h e~^ 2 ) 



4m 2 e -2m 0b v v — ) p"* > 1 2 

There is one trivial solution at b — > oo. In addition, we expect a local maximum and a local 
minimum. These two are at 

v v - v h e- £ ' mob / 2 = 5 (C.7) 

where S is the solution to 

1 + e' 4 1 + e' 2 K J 



1 ' ±:T7 ^Lf 1 + M 1/2 . (c.9) 



That is 



v v 41 + e' 2(1 + e') V 4 
Then at these two extremum 

V"(b ) = -8e'Smlv v e- 2mobo + 0(e /2 ) . (CIO) 

Since this is O(Se'), we keep S to 0{\fd). Using the above solution for 5, we find that to this 
order 



v v . (C.ll) 



So 5 is the solution corresponding to the local minimum. This is expected, since there 
should be local maximum between this local minimum and the minimum at b — > oo. Then 

V"(b) = Ae' 3/2 m 3 v 2 e- 2mobo + 0{t' 2 ) . (C.12) 

Using Eq. ( f4.11|) and VLq = e _mof> °/ 2 , the radion mass is 

2 _ 8 3/2 m v 2 2 

~ e 3/2 A 2 w , (C.13) 

since tuq ~ Mpi and v v ~ Mp{ 2 . So the radion mass is O(TeV) 2 . For the GW solution, 
e' ~ 1/40, so the radion is roughly an order of magnitude below A^. 
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D Newton Force Law Between Particles on the TeV 
and Planck Branes 



We have seen that the expansion in the effective four dimensional theory is given by 

R2 = | = ^?r^| 0* + + pm») ■ (D.i) 

That two mass scales, O(Mpi) in p*, and O(TeV) in pfig, appear in the expansion rate for 
the effective theory is at first surprising. This suggests that observers living on the TeV 
brane will measure "Planck-massed" particles living on the Planck brane to really have 
"Planckian" masses, even though the cutoff on the TeV brane is O(TeV). In fact, this is 
confirmed by a computation of the Newtonian force law between a particle on the Planck 
brane and a particle on the TeV brane. To see this, the 5D Lagrangian for a point mass 
particle with bare mass mi living on the Planck brane, and a point mass particle with bare 
mass m 2 living on the TeV brane is 

-mi J d 5 x5(y)5\x - z x {r))^g~^x~± v - m 2 J d 5 xd(y - l/2)5 4 (x - z 2 {r))^x~± v ,(D.2) 

where g^ v = Q 2 (y)g lll/ and indices are raised and lowered with respect to g. The stress-energy 
tensor (with respect to g) for one of these particles is 

= -mi^—5(y - yi)S 4 (x - z,i{r))—^^= = rmS^g^) . (D.3) 
V9 Jg^XuXv 



But if we express S in terms of g, and raise and lower indices with g, then 

Tfj, v = rriitt^ty^S^ig^). (D.4) 
But with g = Q 2 + h, the (linearized) interaction with gravity is 

J dy^%M v = m l J dytf^gS^M 1 ' ■ (D.5) 

This should be expressed in terms of the graviton zero mode, which is h/j, u (x, y) = h^ u (x) 
where h^ u (x, y) = Q 2 (y)h fiu (x, y). Then the coupling of the particles to the zero mode is 

-rmj dyQJgS^x, y)h^(x) . (D.6) 

So we see that each particle couples to h with strength The Newtonian force 

computed from the exchange of the zero mode only, h(x,y) = h(x), is then 

F N = G N m 1 Vt(y 1 )m 2 Vt(y 2 ) / r . (D.7) 

So the Newtonian force law and Hubble Law for matter living on different branes are con- 
sistent. This result is consistent with [IB , where the authors found that the physical masses 
are given by Q(yi)rrii. Therefore both mass scales appear in the 4D effective theory. 
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E Couplings of Bulk Matter 



In this Appendix, we show what the natural suppression scale for couplings of bulk fields 
are to matter on the Planck or the TeV branes. Following ]T8|, we consider a general bulk 
scalar $ (not necessarily the GW field of radius stabilization) with bulk action 



5 $ = ~ J d b x^G (G AB d A §d B <$> - m 2 & 



(e.i; 



where Gab is the full 5D metric which includes the RS warp factor, and m is the bare 5D 
mass of this scalar, assumed to be of the order of the 5D Planck mass M. We assume that 
this field $ has some generic couplings with matter on a brane of the form 



ind 



M 3 p Mi '' 



;e. 2 ) 



where A+q denotes a composite operator of fields living on one of the branes, and has 
mass dimension 4 + q. We assume that in the bare Lagrangian every suppression scale 
is proportional to the M. Performing the usual field redefinition for the fields living on 
the brane results in O i+q — *> Q~ q (yi)0 4+q , where O i+q has canonically normalized fields 
and physical masses, and Q{yi) is the warp factor at the position of the brane, Q(0) = 1, 
= 1CT 15 . We can see that the effect of this this redefinition is to turn the suppression 
factor M q for the operators on the brane into ~ A^/ for matter on the TeV brane. For 
matter on the Planck brane the suppression factor remains unchanged since no rescaling of 
the fields is needed. 

We decompose the field <3> as 



(E.3) 



where 18 



Vn{y) 



n(y) 



-2 r 



N„ 



\m 



b Y 



m 



m 



(E.4) 



J v and Y v are Bessel functions of order v 



4 + 



m„ is the mass of the Kaluza-Klein 



mode \l/ n , and N n and b nu are normalization constants. With this substitution for $ the 
interaction of a particular KK mode to the brane fields is given by operators of the form 



M 3 p (tt( yi )M) q 



(^(^)) 2P - 



(E.5) 



Thus we need to find the approximate value of the wave function of the KK modes at the 
position of the two branes. First we calculate the suppression at the TeV brane. We will use 
the approximate values v ~ 2, N n ~ Qq 1 . For small KK masses b nv is approximately given 



31 



by 7rm^/4mQ ~ Qq, thus for arguments of order one in the Bessel functions the contribution 
of Y v can be neglected. Thus the approximate value of <p n is given by [] 

Ml/2) ~ ^ = ^ (E-6) 

which then turns every factor of M in ( |E.5| ) into ~ A^. Thus the interactions with matter 
on the TeV wall are suppressed by the TeV scale as expected. 

To find the value of the wave functions close to the origin, we expand ^(^n) ~ tfi 
and Y 2 (x n ) ~ — \ for small values of x n = m n /m . Thus now the contribution of Y 2 will 
dominate, and we get the approximate value 

^n(O) ~ r^X^ ~ Q Q . (E.7) 
"0 X n 

Thus the interactions on the Planck brane are suppressed by a scale higher than Mpi, namely 
MpilO 15 . This can be understood by the suppression of the wave functions of the KK modes 
close to the Planck brane. If there was no suppression of the wave functions, then the 
couplings would be suppressed by Mp[. The small values of the wave functions yield an 
additional suppression, thus resulting in this highly suppressed interaction strength. 

In fact, this wavefunction suppression at the Planck brane is required to maintain the 
large hierarchy at the loop level in the effective 4D theory. Consider for example a scalar 
field living on the Planck wall with mass of O(Mpi), which has renormalizable couplings 
to scalar fields in the bulk. In the 4D effective theory, the bulk scalar fields appear as a 
tower of KK modes, with masses starting at O(TeV), whereas the Planck scalar has O(Mpi) 
mass. Without any wavefunction suppression of the bulk scalars at the Planck wall, loops 
of this Planck scalar field generate 0(M Pl ) quadratic divergences, which would destabilize 
the hierarchy. With a wavefunction suppression of Qo , however, the quadratic divergences 
are O (TeV) 2 and thus do not destabilize the hierarchy. 
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